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Abstract 

The phenomena of giant halo and halo of neutron-rich even-Ca isotopes are investigated and 
compared in the framework of the relativistic continuum Hartree-Bogoliubov (RCHB) and non- 
relativistic Skyrme Hartree-Fock-Bogoliubov (HFB) calculations. With two parameter sets for 
each of the RCHB and the Skyrme HFB calculations, it is found that although halo phenomena 
exist for Ca isotopes near neutron drip line in both calculations, the halo of the Skyrme HFB 
calculations starts at a more neutron-rich nucleus than that of the RCHB calculations, and the 
RCHB calculations have larger neutron root-mean-square (rms) radii systematically in A'^ > 40 
than those of the Skyrme HFB calculations. The former difference comes from difference in shell 
structure. The reasons for the latter can be partly explained by the neutron 3sx/2 orbit, which 
causes more than 50 % of the difference among the four calculations for neutron rms radii at ^^Ca. 

PACS numbers: 21.10.Dr, 21.10.Gv, 21.60.Jz, 27.40.+Z, 27.50.+e 
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I. INTRODUCTION 



The study of exotic nuclei with extreme N/Z ratios has attracted worldwide attention 
since the first discovery of halo in ^^Li Usually for exotic nuclei, their Fermi surfaces are 
very close to the continuum threshold, and the valence nucleons could be easily scattered 
to the continuum states due to the pairing correlations 0, 0]- Thus, theories which can 
properly handle the pairing and continuum states are needed to describe the properties of 
exotic nuclei and to understand the interference of the continuum-energy states to bound 
many-body systems and effects of large neutron excess. 

It is well known that the BCS approximation for the pairing correlations suitable for 
stable nuclei cannot give correct wave functions near the drip lines, see e.g. |^]. Taking 
into account the pairing correlations by the Bogoliubov transformation, relativistic and/or 
non-relativistic mean field approaches have been extensively used to describe and predict 
halo phenomena in exotic nuclei. 

The status of the studies of the halo using the relativistic continuum Hartree-Bogoliubov 
(RCHB) or the non-relativistic Hartrcc-Fock-Bogoliubov (HFB) method include: the pre- 
dictions of halo in the Ne ^J, Na JL], Ca [1] and Zr isotopes near the neutron drip 
line. More details can be found in recent review papers, e.g., 113. As there are more than 
two particles in the weakly bound or the positive-energy region in Ca and Zr isotopes near 
the neutron drip line, the halo phenomena for Ca and Zr are addressed as giant halo fli0|. 
The proton-magic nuclei from O and Pb nuclei are also systematically surveyed, and it is 
found that the halos signature in Zr is the clearest in Ref. There is no evidence for halo 
in Ni, Sn and Pb isotopic chains, similarly as in Ref. T^ |. 



In the non-relativistic HFB approaches, Skyrme calculations with the parameter set SLy4 
predicts the halo in Sn and Ni It is also claimed that the pairing gaps have an effect 
to reduce the halo at the neutron drip line [l^ [13 • 

For comparison with the halos predicted in relativistic approaches [1,1^, ^|, it is necessary 
to investigate the Ca and Zr in non-relativistic approaches, which is very important but still 
missing in the literature. Therefore the present paper is dedicated to investigate the neutron- 
rich even Ca isotopes up to the neutron drip line. The results from both relativistic and 
non-relativistic mean field approaches taking into account the pairing correlations by the 
Bogoliubov transformation will be presented and compared in detail. 

The numerical details used in both calculations are explained in Sec. m and then the 
neutron-number dependence of the energy and the rms radii is examined in Sec. IIIII The 
properties of single particles are discussed in Sec. |TV|in order to understand the halo struc- 
ture. A brief summary is given in Section IVl 



II. METHOD OF CALCULATION 

The detailed formulations of the RCHB method and the HFB method can be found 
respectively in Refs. and [1^, Both calculations are based on the coordinate 

representation and use a box size 20 fm with the spherical symmetry assumed and a mesh 
size 0.1 fm. The quasiparticle states are obtained up to 120 MeV (165-190 radial wave 
functions in ^^Ca), and all of these states are used for calculating potentials. The maximum 
angular momentum of the quasiparticles jmax is y (RCHB) and ^ (Skyrme HFB). The 
parameter sets NL-SH [3 and PKl [19j ( SkM* and SLy4 |2l| ) are used for the RCHB 
( Skyrme HFB ) calculations. For the pairing part, a surface-type delta interaction |^ 1^ 
is used with the density parameter po = 0.152 fm~^ and 0.160 fm~'^ for the RCHB and the 
Skyrme HFB calculations, respectively. 
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FIG. 1: Binding energy per nucleon B/A of even Ca in RCHB calculations with NLSH and PKl 
and the Skyrme HFB calculations with SkM* and SLy4 compared with the experimental data [2^ 
at A'' = 30 and 32. For the details, see the text. 




FIG. 2: Two-neutron separation energy S2n of even Ca in RCHB calculations with NLSH and 
PKl and the Skyrme HFB calculations with SkM* and SLy4 compared with the experimental data 
0. Note that 52„ of SLy4 are very close to those of PKl in 42 < iV < 48. 



The strength of the pairing interaction in mean-field calculations is normally deter- 
mined so as to reproduce the pairing gaps obtained from odd-even mass differences (see 
e.g. Ref. While for nuclei close to the neutron drip line, the pairing energy obtained 

using the Gogny interaction in the pairing channel is adopted to fix the strength of surface- 
type delta pairing interaction due to the missing of experimental data. The Vq is —325 
MeVfm^ in the RCHB calculations, and that of the Skyrme HFB calculations is —365 
(—330) MeVfm^ in the SLy4 (SkM*) calculations. The average pairing evergy of the PKl 
and NL-SH calculaitons —8.2 MeV at ^^Ca is used for determining the Vq of the Skyrme-HFB 
calculaitons. We checked that when jmax = y was used in the Skyrme HFB calculations, 
the rms radii did not change with the re- adjusted Vq from the pairing energy. 



III. NEUTRON-NUMBER DEPENDENCE OF THE HALO STRUCTURE 

Figures H and m show the binding energy per nucleon B/A and the two-neutron separation 
energy in RCHB calculations with NLSH and PKl and the Skyrme HFB calculations 
with SkM* and SLy4, respectively. The most important differences between these calcu- 



3 



4.8 



4.6 



4.4 



4.2 



3.8 



3.6 



3.4 



3.2 



NL-SH 

PK1 ■■■■ 

■ SkM* □ 

SLy4 A... 



Ca 



/T 

// / 

■■'.t' El' 



neutron 



proton 



30 



40 



50 



60 



FIG. 3: The neutron and proton rms radii f„ and fp of even-Ca isotopes in RCHB calculations 
with NLSH and PKl and the Skyrme HFB calculations with SkM* and SLy4 in the neutron-rich 
region. 



lations are location of the two-neutron drip line and a shell gap at iV = 40. The RCHB 
calculations show the two-neutron drip line at = 50 (PKl) or 52 (NL-SH), and the pre- 
diction of the Skyrme HFB calculations is = 46 (SLy4) or 58 (SkM*). The RCHB and the 
SLy4 calculations have clear the shell gap at = 40 (Fig. |21), while the SkM* calculation 
does not show at all. The SkM* calculation has a small extra downward behavior of S2n at 

= 50, however, it is too small to call a shell gap. Therefore, no new shell gap is predicted 
in A^ > 40, in any of the calculations. The two RCHB calculations give the results rather 
close to each other both in B/A and S2n, and the two Skyrme HFB calculations have a 
difference of 0.1-0.2 MeV in B/A. The S2n of SLy4 is close to those of the SkM* calculation 
in 34 < A^ < 40 and then coincides with that of PKl in 42 < A^ < 48. The experimental 
data are available up to A^ = 32 currently [2^. The RCHB calculations reproduce well 
the data of B/A at A^ = 30 (8.550 MeV) and 32 (8.396±0.013 MeV) and S2n at N = 32 
(9.081±0.699 MeV). For S2n a.t N = 30, the SLy4 calculation gives a value very close to the 
measured one (11.499±0.008 MeV). 

Figure 01 shows systematics of the rms radii of neutrons f„ and protons fp of even-Ca from 
A^ = 30 to the two-neutron drip lines. The increase in the curvature of indicates the halo. 
(For attempts to define the criterion of halo or the halo size, see Refs. 05^?^] and references 
therein.) For Ca, the number of nucleons in the positive-energy region of ^^^''^Ca is 0.6- 
2.2, of which average is 1.7, in the NL-SH calculation Js!], and the corresponding value of the 
SkM* calculation is always smaller than 0.5 (the average 0.27). As there are more than two 
nucleons in the weakly bound orbits and the positive-energy region for these nuclei and also 
halo phenomena in neighboring nuclei with incremental neutrons, the halo in Ca isotopes 
are refereed as giant halo, as in Refs. P,l9j. The SkM* calculation has the halo from A^ = 52, 
and the RCHB calculations show the gradual occurrence of the giant halo. Apparently, the 
starting nucleus of the halo of SkM* corresponds to the extra lowering of the S2n- The SLy4 
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FIG. 4: The difference between neutron and proton rms radii — fp in RCHB calculations with 
NLSH and PKl and the Skyrme HFB calculations with SkM* and SLy4 for even-Ca isotopes with 
N > 40. 
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FIG. 5: The difference of neutron rms radii between neighboring even-Ca isotopes with > 40 
fn{N) — fn{N — 2) in RCHB calculations with NLSH and PKl and the Skyrme HFB calculations 
with SkM* and SLy4 in the neutron-rich region. The dashed line was obtained from f„(A^) = 



calculation does not have the halo, because the particle-stable region ends at iV = 46. The 
large difference between f„ and fp shown in Fig. El is identified with neutron skin in the 
region with no halo. f„ — fp is displayed in Fig. |3J in which the differences between the 
parameter sets are clear; the curve of SLy4 is parallel to that of SkM*, and PKl has a larger 
curvature than the other parameter sets in < 50. An important aspect of the halo is 
how f„ changes when two more neutrons are added, thus it is worth showing the difference 
fn{N) — fn{N — 2) (Fig. El). When there is no halo, two additional neutrons increase f„ by 
0.05 fm, and the increase rises to 0.11 fm in the halo nuclei. In all of Figs.OHSl the difference 
in fn is clear between the RCHB and the Skyrme HFB calculations. It is to be noted that 
fp is not constant, reflecting the self-consistency of the neutrons and protons interaction. 



IV. SINGLE-PARTICLE STRUCTURE 

To understand the halo phenomena and the different predictions in RCHB and Skyrme 
HFB calculations, the nucleus ^^Ca is taken as an example. Figure El depicts the single 
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FIG. 6: The neutron sing le particle level e^f" for «'^Ca in RCHB calculations with NLSH and PKl 
and the Skyrme HFB calculations with SkM* and SLy4. 

particle level e'^'^ in RCHB calculations with NLSH and PKl and the Skyrme HFB calcula- 
tions with SkM* and SLy4 for neutrons in ^^Ca, i.e., the diagonal elements of the mean- field 
Hamiltonian in the canonical basis. 

The order of the levels is essentially the same for all the calculation. Furthermore the 
spectra for RCHB with PKl and NL-SH are quite similar. However, the level lsi/2 in SkM* 
is apparently higher than those of the other parameter sets. From left to right, the shell gap 
become less apparent. The shell gaps at = 2, 8, 20, 28, and 40 are clearly shown in the 
RCHB spectra, while the shell gaps of = 28 and 40 do not appear in that of the SkM* 
calculation (see also Fig. |2I). 

It has been pointed out that the neutron 3si/2 level played an important role in the 
giant halo. Thus, comparison is made in Fig. [7| for the rms radii of each single-particle 
orbit (canonical basis) in RCHB calculations with NLSH and PKl and the Skyrme HFB 
calculations with SkM* and SLy4. The lines of versus e'^^^ are quite similar to each other 
except for the lsi/2 and the region around e'^^'^ = 0. This indicates the similarity of the 
potential wall except for the bottom of the potential wall. A magnification around e^'^^ = 
is displayed in Fig. [7|d. The most significant difference is 3si/2: ^3si/2 of PKl is more than 
2 fm larger than that of SkM*. These differences are clearer in Fig. |8l which illustrates as 
radii of circles. 

Table m shows e"^"^, f^, occupation probability f^, and ^v'^fa for the neutron 3si/2 orbit, 
as well as f„ of ^^Ca. Both 7^3^1/2 and "^3^1/2 ^^ake f„ of PKl the largest one. Since the 
contribution of the orbit to f„ is given by jfVlgi/2^3si/2, it is seen from Tab. H] that 54 % 
(0.082 fm) of the difference f„(PKl) — fn(SkM*) = 0.153 fm comes from the neutron 3si/2. It 
is also noted that the higher energy does not necessarily mean the larger spatial distribution. 

The difference in the neutron 3s 1/2 orbit is shown also in Fig. IHl which illustrates the 
radial wave functions squared. The tail of PKl has appreciably longer distribution than the 
others, and as a counter part, the amplitude of PKl is smaller in the inner region. The 
number and locations of the nodes are the same for the four lines. The neutron density 
distributions have appreciable difference (Fig. [TUI) in accordance with the behavior of Fig. El 

Figure im depicts central potentials. These are the summation of the scalar and the vector 
potentials in the RCHB calculation, and the summation of the to and terms in the particle- 
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FIG. 7: a) The rms radii of each neutron single-particle orbit in canonical basis Tq, as functions 
of its single-particle energy e^^'^ . b) Same as a) but for the region around e'^^'^ = 0. 



TABLE I: Properties (ej^f°, Tq,, occupation probability u^, and 
jj'v'^fa ) of the neutron 3si/2 orbit and Vn of ^^Ca. 





^can 
^3sl/2 


^^351/2 


^3sl/2 
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[MeV] 


[fm] 




[fm] 


[fm] 


PKl 


0.133 


8.144 


0.239 


0.0846 


4.363 


NL-SH 


0.640 


7.240 


0.089 


0.0280 


4.314 


SLy4 


L939 


6.856 


0.016 


0.0043 


4.205 


SkM* 


1.445 


5.765 


0.010 


0.0025 


4.210 



hole potential in the Skyrme HFB calculation, see Eq. (A. 5a) of Ref. [1h]. A significant 
difference is that the SkM* calculation has smaller depth than the other calculations, as is 
anticipated from Fig. and another is a slight difference in the tail region (see the inset of 
Fig.lll|). It should be noted that the tails of the RCHB calculations have longer distributions 
than the Skyrme HFB calculations, and PKl has the longest one. Figure ^1 displays the 
pairing potential 

where p(r) is the total nuclear density, and Pnir) is the pairing density jl^ of neutron. 
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FIG. 8: The rms radii of each neutron single-particle orbit in canonical basis in RCHB calcu- 
lations with NLSH and PKl and the Skyrme HFB calculations with SkM* and SLy4 displayed as 
radii of circles. The proton one fp (the blue center circle in the on-line version) is included for a 
comparison. The grade of color corresponds to the order of f^. 



For the RCHB calculations, the S = 1 component of the pairing interaction is removed 
explicitly. Again, the line of PKl has the longest tail, therefore, the large f3si/2 of PKl can 
be understood in terms of the tails of the central and the pairing potentials. 

Finally, we mention the importance of the pairing correlations. It turned out that when 
the strength of the surface-type pairing interaction was determined so as to reproduce the 
neutron pairing gaps of ^^'^^'^^'^°Ca, 1.0-1.9 MeV, obtained from the odd-even mass differ- 
ences of the experimental data the pairing energy of the SkM* calculation was —(35-39) 
MeV in 40 < < 50, and the kink of f„ at = 50 was completely smeared out (the circles 
in Fig. IT^ . It is noted that the pairing energy is more than 3 times larger than the values 
of the RCHB calculations 8-9 MeV. On the other hand, the volume-type pairing interaction 
Vo5(r" — r') can reproduce the pairing energy of the RCHB calculation and the experimental 
pairing gaps simultaneously, and the f„ has the kink (the triangles in Fig. 1131) . Therefore, 
the kink is strongly influenced by the pairing interaction. (See also Refs. ^13«. il4|-) 



V. CONCLUSION 

In this paper, the phenomena of the giant halo and halo of the neutron-rich even-Ca 
isotopes have been investigated and compared in the framework of the RCHB and the 
Skyrme HFB calculations. With the two parameter sets for each of the RCHB and the 
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FIG. 10: Neutron density distribution of ^^Ca calculated with the four parameter sets. 

Skyrme HFB calculations, it has been found that although the halo phenomena existed for 
Ca isotopes near the neutron drip line in both calculations, the halo of the Skyrme HFB 
calculations started at a more neutron-rich nucleus than that of the RCHB calculations, and 
the RCHB calculations had larger neutron rms radii systematically in A'" > 40 than those 
of the Skyrme HFB calculations. The former difference comes from difference in the shell 
structure. The reasons for the latter can be partly explained by the neutron 3si/2 orbit, 
which causes 50 % of the difference in the neutron rms radii among the four calculations at 
66Ca. 
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FIG. 11: Central potential of the neutrons of ^^Ca. Inserted is a magnification of the surface 
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FIG. 12: Pairing potential of the neutrons of ^^Ca. 
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FIG. 13: f„ of the SkM* calculations with three pairing interactions: the surface-type pairing 
with the strength determined from the pairing energy of the RCHB calculations (squares), the 
surface-type with the strength determined from the pairing gaps (circles), and the volume-type 
determined from the pairing gaps (triangles). The lines of the second and the third ones are very 
close to each other. 
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